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Abstract 

The most common decay used for measuring 2/3 s , the phase of B®-B® mixing, is 
B® — > J/ip(f). This decay is dominated by the colour-suppressed tree diagram, but 
there are other contributions due to gluonic and electroweak penguin diagrams. 
These are often referred to as "penguin pollution" (PP) because their inclusion in 
the amplitude leads to a theoretical error in the extraction of 2/3 s from the data. In 
the standard model (SM), it is estimated that the PP is negligible, but there is some 
uncertainty as to its exact size. Now, <f> c s cs (the measured value of 2/3 s ) is small, in 
agreement with the SM, but still has significant experimental errors. When these 
are reduced, if one hopes to be able to see clear evidence of new physics (NP), it is 
crucial to have the theoretical error under control. In this paper, we show that, using 
a modification of the angular analysis currently used to measure (f> c s cs in B° s — >■ J/i^cj), 
one can reduce the theoretical error due to PP. Theoretical input is still required, 
but it is much more modest than entirely neglecting the PP. If <p c s cs differs from the 
SM prediction, this points to NP in the mixing. There is also enough information to 
test for NP in the decay. This method can be applied to all B®/B® — > V1V2 decays. 
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In the presence of B®-B® mixing, the mass eigenstates B S L H [L (H) corresponds 
to "light" ("heavy")] are admixtures of the flavour eigenstates B° s and B Q S : 

\Bl, H )=p\B s )±q\B s ) , (1) 

where the complex coefficients p and q satisfy \p\ 2 + \q\ 2 = 1. States which are B° s 
or Bg at t = then evolve in time into an admixture of both states, leading to 
the time- dependent states B®(t) and B^(t). If both B® and B° s can decay to the 
final state /, there is an indirect (mixing-induced) CP-violating asymmetry (CPA) 
between the rates \B° s {t) — > f\ 2 and \B^(t) — > f\ 2 . The indirect CPA measures 

qM 
V A 



Im(^), (2) 



where A{ and A{ are the amplitudes for B® — > f and B® — > /, respectively. A{ can 
be obtained from A{ by changing the signs of the weak phases. 

B° s -B° s mixing is dominated by the box diagram with an internal t quark, so 
that q/p = (V t * b V ts /V t bV t * s ) = exp(2i arg(V^T4 s )). This is phase-convention de- 
pendent. Suppose now that A{ is dominated by a single decay amplitude with 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements V* b V cs . We then have 
A{/A{ = (V cb V* s /V* b V cs ) = exp(2iarg(V r c6 V A c *)), which is also phase-convention de- 
pendent. However, the product of these two quantities is 

q A{ v t tv ts v cb v; s 
paI v; b v cs v tb v t * s 

where 



e 2i/3s , (3) 



Pa = arg 



v; b v cs 



(4) 



This is phase-convention independent, and hence physical. The indirect CPA mea- 
sures sin2/3 s . 

In the standard model (SM), 2(3 S is expected to be very small: the SM prediction 
is 2/3 s = 0.03636±0. 00170 rad pp. For this reason there is great interest in measuring 
this quantity. A value for 2f3 s significantly different from zero would be a smoking- 
gun signal of new physics (NP). 

The main decay mode used for measuring 2(3 S is B° s — > J/ip<j), which is the ana- 
logue of the "golden" mode B^ — > J j ipK s used to measure the B®-B® mixing phase 
2/3. Since B° s — > J/i(j<f) involves two final-state vector mesons, a time-dependent 
angular analysis of this decay (with J/ip — > £ + £~ and <fi — > K + K~) must be used 
to disentangle the CP + and — final states. The most precise measurement of the 
indirect CPA has been performed by the LHCb Collaboration [2]. They finc0 

(f)f s = -0.001 ± 0.101 (stat) ± 0.027 (syst) rad , (5) 



4 j3 s is defined in Eq. (j4]). For the measured B a s -B Q S mixing phase, it is common to use the 
symbol cf) c s cs (or occasionally just 4> s ), which is equal to — 2f3 s in the SM. 
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in agreement with the SM. (Still, the errors are large enough that NP cannot be 
excluded.) 

Now, the above formalism holds for the case where A{ is dominated by a single 
decay amplitude. However, if A{ contains a second decay amplitude with a differ- 
ent weak phase, then the indirect CPA no longer cleanly measures 2/3 s - there is 
a theoretical uncertainty due to the presence of the second amplitude. It is often 
the case that the first and second amplitudes correspond to tree and penguin dia- 
grams, respectively. The theoretical uncertainty is thus generally called "penguin 
pollution." 

For B° s — > J/ip4>, the decay amplitude can be written 

( n /\ 

where X q = V* b V qq i. Here, C , P and P' EW are the colour-suppressed tree, gluonic 
penguin and electroweak penguin diagrams, respectively. (As this is a b — > s transi- 
tion, the diagrams are written with primes.) The index 'i' of P[ indicates the flavour 
of the quark in the loop. The unitarity of the CKM matrix (Au + Ac + A^ = 0) 
can be used to eliminate the t-quark contribution: 

A m* = X (s)(c> + p' ct - 2 -p' ew ) + \^\p' ut - 2 -p' ew ) 

= e ia ^ v '» Vcs) [Ai + e'' 7 ^] , (7) 

where P' ct = P' c —P' t , P' ut = P' u — P[. In the second line, A\ and A 2 are combinations of 
diagrams, and include the magnitudes of the corresponding CKM matrix elements. 
We have explicitly written the weak-phase dependence of the terms in the ampli- 
tude^]. (The phase information in the CKM matrix is conventionally parametrized 
in terms of the unitarity triangle, in which the interior (CP-violating) angles are 
known as a, (3 and 7 [3].) A\ and A 2 contain strong phases. 
We now have 

q _Af^ e2i ,A 1 + e-*A 2 
pAl A l + e^A 2 

The point here is that the A 2 term represents the penguin pollution (PP). There are 
several ingredients in estimating its size. First, |V^* 6 \4 S |/| V^V^I ~ 0.02. Second, in 
Ref. [1], it is estimated that \P' ct /C'\ ~ \Pu t /C'\ = 0(1). However, here the penguin 
diagrams are both OZI-suppressed, so these ratios are much smaller. Third, it is 
expected that \P'ew/C\ = 0(A), where A ~ 0.2 [4]. Putting everything together, it is 
estimated theoretically that IA2/A1J = O(10~ 3 ) [5]. Thus, to a good approximation, 

5 The weak phase in the coefficient of A2 is a,vg{V* b V us /V* b V cs ), and so should technically be 
called 7 S . However, it differs very little from 7 = ^g(V* b V u d/V* b V c d), and so we refer to it by this 
name. 
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the A2 term, i.e. the PP, is negligible. In this case, the indirect CPA in B® — > J/i/j(p 
cleanly measures 2(3 S . 

Indeed, in the determination of <p c s cs [Eq. (jSJ)], the LHCb Collaboration neglected 
the PP [2] (A{ = A{ was assumed). However, the result reintroduces the question 
of the size of the PP. To be specific, we now know that <p c s cs is small. Suppose 
a future measurement finds <j) c s cs = 0.1 ± 0.01 rad (error is statistical only). This 
disagrees with the SM prediction and therefore suggests NP. But perhaps the size of 
the PP has been underestimated, so that there is really a nonzero theoretical error 
associated with this measurement. It is possible that, when this error is taken into 
account, the discrepancy with the SM disappears, so that there is no signal of NP. 
This situation is not at all improbable - the theoretical prediction of O(10 -3 ) for 
the size of the PP includes estimates of the hadronic matrix elements. But these 
are notoriously difficult to determine with certainty. It is not impossible that the 
true size of the PP is larger than its theoretical estimate, and that this introduces 
a theoretical error which is important when the measured value of <fi c s cs is small. 

This issue was raised several years ago in Ref . [B] • There it was suggested that the 
PP term can be measured in the time-dependent angular distribution of the flavour- 
specific b ->• d decay B° s J/ipK*°(-^ n + K~), and then related to B° s J/ip4> 
using flavour SU(3) symmetry. The difficulty here is that the value of SU(3) breaking 
is unknown, and this is extremely problematic given that we need a precise value of 
the size of the PP. 

The purpose of this paper is to point out that, in fact, using the present time- 
dependent angular analysis, the theoretical error due to PP associated with <j) c s cs can 
be reduced. In this method, one retains the PP term. One piece of theoretical input 
is still needed to extract <j) c s cs , but the theoretical error is quite a bit smaller than 
when the PP term is neglected from the beginning. In addition, if it is concluded 
that NP is present, there is enough measured information to determine whether the 
NP is in the mixing and/or the decay. The method is described below (here we 
follow the procedure as described in Ref. [7]). 

Given that B° s — > J/ip<f) has two final-state vector mesons, its amplitude can 
be separated into 3 helicities: the polarizations of the vector mesons are either 
longitudinal (Ao), or transverse to their directions of motion and parallel (An) or 
perpendicular (A±) to one another. In addition, since the <fi is detected through its 
decay to K + K~, one must also allow for the possibility that the observed K + K~ pair 
has relative angular momentum / = (S'-wave) [8]. That is, the angular distribution 
of the decay B® -> J/if>(-> £ + £~)<t>(-> K + K~) involves the 4 helicities 0, ||, _L, S. 
In going from the decay to the CP-conjugate decay, Ah —> r/hAh (h = 0, ||,_L, S), 
in which the Ah are equal to the Ah, but with weak phases of opposite sign, and 
Vo = V\\ = +1, V± = ri ls = -1- 

The angular distribution is given in terms of the vector Q = (8, ijj, <p) in the 
transversity basis [HI EDI ED]- This basis is defined as follows: in the J/i/j rest frame, 
one has a righthanded coordinate system in which the x axis is parallel to and 
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the z axis is parallel to px- x Pk+- In this frame, 6 and tp are the azimuthal and 
polar angles, respectively, of the The angle ip is the angle between p K - and p j/^ 
in the rest frame of the <fi. We have 

— 1 ' oc £ fc fc (t)/ fc (n) . (9) 

dtd\ I Tr x 

The explicit expressions for the are given in Ref. [7]. 

The hk(t) are determined as follows. The time dependence of the helicity ampli- 
tudes Ah (h = 0, ||, _L, S) is given by [12] 



A h (t) = g + (t)A h (t = 0) + Vh Z g.(t)A h (t = 0) . (10) 

p 

g±(t) contain the information about the time evolution of B®(t) and B^(t). We use 
the following established relations involving g±(t): 

\g±(t)\ 2 = ie" rs ^cosh(Ar s /2)t ±cosAm s t) , 

g* + (t)g-(t) = \z~ Vt (y ~ sinh (AT s /2)t + isinAm s t) . (11) 

The hkit) are given by 

h(t) = \A (t)\\ h 2 {t) = \A^t)\\ h 3 (t) = \A x (t)\\ h 7 {t) = \A s {t)\\ 

h 4 (t) = Im (A ± (t)Al(t)) , h G (t) = lm(A ± (t)A* (t)) , 

h 5 (t) = Re(A (t)Al(t)) , h 9 (t) = lm(A ± (t)A* s (t)) , 

hs(t)=Re(A\\{t)A* s (t)) , h 10 (t)=Re(A {t)A* s (t)) , (12) 
and can be written as 

+ a k cosh (Ar s /2)t + b k sinh (Ar s /2)t] . (13) 

By measuring the time-dependent angular distribution and fitting to the four time- 
dependent functions, T s and Ar s can be determined, as well as the coefficients a k -d k . 
(Note that the width difference in the B®-B® system, AT q , is sizeable only for B° s 
decays. Thus, this method for dealing with PP cannot be used for B® decays.) 

There are also the functions hk{t) in the angular distribution of the CP-conjugate 
decay B® — > J/xfxf). They can be obtained from the hk(t) by making the following 
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changes: Ah -H- rjhAh and 4> c s cs — > —4>f s . Each of the a^-d^ of the hk(t) is the same 
as the dk-dk of the /ifc(t), up to a possible sign. Both hk(t) and /ifc(t) should be 
included in the fit. 

Also, the hk(t) contain both CP-conserving and CP- violating terms. However, 
when one calculates the sum or difference of hk(t) and hk(t), certain terms cancel, 
so that these sums and differences are purely CP-conserving or CP-violating. In 
particular, the untagged sums h k (t) + hkif) are CP-violating for k = 4,6,8,10. 
These correspond to triple-product terms. Similarly, the differences hk(t) — hk{t) 
(k = 1, 2, 3, 5, 7, 9) are CP-violating. These can only be constructed through tagged 
decays, and contain the direct and indirect CPA's. 

Now, the ultimate goal is to measure <^ cs . To this end, it is necessary to express 
the cik-dk in terms of all the theoretical unknowns in order to determine what exactly 
can be extracted. The bk and dk involve terms of the form 

Re(^A* h A h ^\ , lm(U* h A h ^ . (14) 

The helicity amplitudes can be multiplied by an arbitrary phase, so that the form 
of the individual terms is uncertain. In order to fix the phases, in what follows we 
adopt the convention that the dominant SM decay amplitude has no weak phase. 
That is, the phase term e~ 2ia,rg ( v * bVcsS> is factored out and combined with the phase 
of q/p to make 2f3 s , as in Eq. ([8]). Thus, any terms of the form Ke(q/p) or lm(q/p) 
depend only on the mixing phase f3 s . 

In general, the a^-d^ are expressed in terms of 16 unknown parameters: the 
magnitudes of the A h and A h (8), their relative phases (7), and 0^ cs . For the phase 
differences we define 



Sij = arg(Aj) - aig(Aj) , 
5ij = arg(Ai) - arg(Aj) , 

Dij = arg(A) - arg(Aj) , (15) 

where i,j are any of the 4 helicities 0, ||, J_, S. Of course, there are many phase 
differences, but only 7 are independent. It is straightforward to show that 



-Dij Dji -Da ~~\~ , 

D u - Djj = 5ij - . (16) 

Suppose first that only the dominant SM amplitude is retained (the "1-amplitude 
method") - this is A\ of Eq. ([7]). In this case we have Ah = Ah (h = 0, ||,_L, S), 
so that Sij = = S^. Now there are 8 unknown parameters - the magnitudes of 
the Ah (4), the relative strong phases (3), and <j) c s cs . (In their analysis, LHCb fix the 
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normalization \Aq\ 2 + |A||| 2 + |A^| 2 = 1-) The coefficients ak-dk are given in Table 
[TJ It is clear from this Table that all unknown parameters can be determined from 
the measurements of ak-dk- In addition, there is a great deal of redundancy, which 
allows for a reasonably precise determination of these unknowns. Indeed, this is the 
method used by the LHCb Collaboration [21 [7J. However, as noted above, there is 
a potential theoretical error associated with the neglect of the PP. 



ak 


hi 




2\A 


h 2 




2\A\\ 


h 3 




2\A ± 


hi 









2|A| 




h 6 







h 7 




2|A 5 


h s 







/>!) 


2\A ± 


Us\ 


hw 








h Ck dk 

P -2|A | 2 cosc/>f s 2\ A \' 2 sin 4>f s 

| 2 -2|il||| 2 cos0f 2|A||| 2 sm<^f s 

| 2 2\A ± \ 2 cos (j)f s -2|A_L| 2 sin(#p 

-2| J 4j_||^4|||cos(5_ L || sin^f 5 2|A_|_||A||| sin^y -2|AjJ|A||| cos<5_l|| cos</>f s 

cos J|| — 2 1 An ||A | cos Sun cos <^ cs 2|Aii ||A | cos 5\\ sin<^ cs 

-2\A ± \\A \cos8± sm(j)f s 2|A_|_||A | sin<5_Lo -2\A ± \\A \cosS± cos<f)f s 

| 2 2|A s | 2 cos^f s -2\A S \ 2 sin <f>f 

-2| Ay 1 1 A s | sin 6\\ S sin 0f 2| Ay 1 1 A s | cos -2| Ay 1 1 A s \ sin 5|| S cos <pf s 
sm5± s 2|A_l||A s | sin <5 ±iS cos 0f s -2|A_l||A 5 | sin 5± s sm4>f s 
-2|A ||A s |sm(5ossm^f s 2|A ||A g | cos 5 0S -2|A ||A g | sin 5 0S cos(j)f s 



Table 1: Coefficients ak-dk of Eq. ( {IB"]) in terms of \Ah\ (h = 0, ||,JL, S) and <5y 
[Eq. ( fl5l) ] for the case where A/i = A^. 

Suppose instead that no assumptions about the decay amplitude are made. That 
is, we allow for the possibility of more than one contribution to the amplitude (the 
"2-amplitude method"). This can arise through the inclusion of the PP term and/or 
a NP amplitude. The upshot is that Ah ^ Ah, which is the most general possibility. 
Here the coefficients a^dk of Table [I] are modified - in particular, the and 
which vanish in Tabled] are now nonzero. The a^dk for this case are given in Tables 
[2]|D All of these quantities can be measured in the time-dependent angular analysis. 

Now there are 16 unknown parameters: the magnitudes of the Ah and Ah (8), 
their relative phases (7), and (f> c s cs . \Ah\ and \Ah\ can be obtained from ak and Ck for 
k — 1,2,3, 7 (Table [2]). The 3 independent relative phases among the A h (the Sij) 
and the 3 independent relative phases among the Ah (the dy) can be found from 
and Ck for k = 4-6 and 8-10. For the missing relative phase - a Dij - one has to look 
at the bk and dk (Tables [3|4|) . However, all terms in these tables are a function of 
(jfcs _ Thus, it is not possible to separate 0^ cs and the solely from the data. 
Even in this case, where the PP has been retained, theoretical input is necessary. 
(Note that this not an accident - it is a direct consequence of the fact that one 
requires a convention to fix the phases in Eq. (1141) .) 

On the other hand, it is obvious how to formulate the theoretical input - one has 
to choose a value for one of the D^. For example, suppose we set D 00 = 0. With 
this input, we can now obtain (f> c s cs from a fit to the data. Note that no assumptions 
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(gfc + c fe )/2 (o fe - c fc )/2 

/n l^ol 2 \M 2 

ho U||| 2 U||| 2 



^3 \A ± \ 2 _ \A_ ± \ 2 

hi \A_i\\A\\\ sinful — | v4j_ | \A\\ | sin <5_i_m 

/z 5 |A||||A | cos (5|| |^4||||v4 | cos (5|| 

h$ \A±\\Aq\ sin5j_o — \A±\ \Aq\ sin S± 

h 7 \A S \ 2 _ \A S \ 2 

h$ \A\\\\ A s\ cos S\\s —\A\\\\As\ cos 5\\s 

h 9 \A ± \\A S \ sm5 ±s \A±\\As\ sm5± s 

h w \A \\A S \ cos5qs -\A \\A S \ cos^ 

Table 2: Coefficients and of Eq. (Tl3l) in terms of and \ Ah\ {h = 0, ||, _L, S*), 
and 5y and [Eq. (TT5|) ]. 



are made about the magnitudes and other relative phases of the A^ and A^ - they 
are also obtained from the fit. 

Now, the theoretical error on (p c s cs is equal to the theoretical error associated with 
the assumption that Dqq = 0. This can be computed in the SM. We write 

A = A 1>0 + ^A 2 , , (17) 
where the A 2t o term represents the PP. Then 

Ax) = arg(A*A ) ~ 2 ^4 cos A sin 7 , (18) 

|Ai,o| 

where A is the strong-phase difference between A 2 ^ and A lj0 . We therefore see that 
the theoretical error associated with D 00 = is 0(|A 2j o|/|^4i,o|)- We can get a bit 
more information from experimental measurements: 



2 I ^2,0 1 

\Ai,o\ 

The theoretical error on df™ can then be written 



|4> 


2 _ 


A 


2 r 


\Ao 


2 + 


A 


2 ai 



sin A sin 7 . (19) 



A0f s ~ cot A — . (20) 

Thus, to get a numerical value for A<p c s cs , we need a theoretical prediction for the 
strong-phase difference A. 

We therefore see that the 2-amplitude method requires theory input, with an 
associated theoretical error of 0(|A2,o|/|-4i,o|)- At first glance, this does not seem to 
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h 

h -2|A ||A |cos(</>f -D 00 ) 

h 2 -2\A\\\\_A\\\ cos(0f s - D|m) 

ha 2|A ± ||A ± |cos(^f -D ±± ) 

K sin(0f s - - |4±||A||| sm((pf s - D M \) 

h 5 -\A {l \\A \ cos(0f s - Do,,) " l^llll^ol cos(0f s - D|, ) 

/i 6 -|A ± p | sin(0f - D 0± ) - |AJ|A)| sin(0f - D ±0 ) 

/i 7 2|A s ||^|cos(^f 



/>8 


-l^lll 


\M 


cos( 




-D\\s) 


+ 




1^1 


cos(i 




-D S \\) 


/lg 




\M 


sin(c 




-D ±s ) 


+ 




\A S 


sin( 








-\M 


\M 


cos( 


0f 


-D os ) 


+ 


l^ol 


\M 


cos(i 




-Dso) 



Table 3: Coefficients bk of Eq. (fT3|) in terms of \Ah\ and (/i = 0, l^-L,^), Dij 
[Eq. ([T5]) ], and 0f s . 



be much of an improvement over the 1-amplitude method. However, let us examine 
the theory errors of both methods in more detail. While the 2-amplitude method 
requires one assumption - D 00 = - the 1-amplitude method in fact requires 8 
assumptions: 

\A h \ = \A h \ (h = 0,\\,±,S) , ^ = 4, (ij= \\0,±0,0S) , ^00 = 0. (21) 

Furthermore, there is a theoretical error associated with each assumption. This is 
due to the fact that, if the true values of the parameters do not obey the assumption, 
the imposition of the assumption will lead to a theoretical error in the extraction 
of <p c s cs . It is true that the theoretical error associated with each assumption is 
only 0(|A 2 /Al|). However, when one takes into account the fact that multiple 
assumptions are involved, the net theoretical error is, in reality, quite a bit larger. 
The theoretical error of the 2-amplitude method could well be an order of magnitude 
smaller than that of the 1-amplitude method. 

In fact, it is possible to quantify this difference by pushing the analysis a bit 
further. Starting with the 1-amplitude method, we can relax an assumption by 
adding a new unknown parameter. For example, suppose we allow \Aq\ ^ \Aq\, but 
the other assumptions of Eq. (I2T1) are retained. Now, when doing the fit to the data, 
there are 9 unknown parameters - the \Ah\ and \Aq\ (5), the relative strong phases 
(3), and <fi c s cs . In addition, the theoretical expressions in Tabled] are modified; the 
correct expressions can be found from Tables [2JIU but imposing the assumptions that 
have been retained. The fit will give the preferred values of \Aq\ and | 1 , so that 
we can see to what extent (if any) the assumption that \Aq\ = \A \ is violated. The 
fit will also give the preferred value of <j) c s cs . This can be compared with the value 
found in the full 1-amplitude method, so this will give the theory error associated 
with this assumption. To be conservative, in computing the theory error, one should 
consider the full allowed (la) range for (f> c ^ s . 
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dk 



h 


2\A \\A \ 


sin(0f - D 00 ) 




2|A||||A||| 


sin(0f s - D|m) 


h 


-2|Aj.||Aj 


J sin(0f - D ±± ) 



h A -|A ± p|||cos(0f s -D\\ ± ) - |AJ|A|||cos(>f s --Dj.||) 

h |A||p o |sin(0f - Don) + miMMtf 8 -Dp) 

h 6 -\A ± \\A \ cos(0f - D 0± ) - \A ± \\Ao\ cos(0f - D> ±0 ) 

fc 7 -2|^ s p s |sin(0f -D ss ) 



h 8 \A ll \\A s \sixy((f>f s - D ]ls ) - |A||p 5 |sin(0f - D S \\) 
h 9 -\A ± \\A S \ cos(0f s - D ±s ) + \A ± \\A S \ cos(0f - D s± ) 
h w \A \\A S \ sin(0f - D os ) -\A \\A S \ sin(0f - D so ) 



Table 4: Coefficients dj, of Eq. (fl3|) in terms of \A^\ and \A^\ (h = 0, ||,JL, 5"), 
[Eq. (USD], and (f)f s . 



Now, LHCb has performed an analogous analysis [21 [33] • They allow for \A^\ ^ 
\Ah\, but take this difference to be helicity-independent. That is, they assume that 

144 = A Vh, (22) 
\A h \ 

and add the unknown parameter A to the fit. This non-equality has the effect of 
making nonzero, but helicity-independent, the at and Ck which vanish in Table [TJ In 
the fit, the idea is to see how much A deviates from 1. They find that this deviation 
can be at most ±5%, and leads to a deviation in (p c s cs of ±0.02 rad. Since this theory 
error is small, one might be tempted to conclude that the original assumption is 
actually rather good. However, one has to be careful here. First, we expect the 
theory error associated with each assumption to be reasonably small, Od^/Axl). 
Second, if ±0.02 rad is roughly the size of the theory error associated with each 
assumption, then the total theory error can naively be as large as ±0.16 rad (8 
assumptions), which is non-negligible. It is therefore preferable to make as few 
assumptions as possible. 

This type of analysis can be performed for all the assumptions, except for Dqq = 
0. This assumption, which remains even in the 2-amplitude case, cannot be tested 
experimentally. In the presence of a nonzero D 00 , what is extracted is not 4>g CS , 
but <j) c s cs — D 00 . However, one can determine the theory error associated with all 
the other assumptions. On the other hand, it is clear that it is best to relax as 
many assumptions as possible. That is, rather than assume a particular value for an 
unknown parameter, one adds this parameter to the fit, and uses the best-fit value. 
There is still an uncertainty in <p c s cs associated with the allowed ranges of these 
new parameters, but this is statistical, and not systematic. Ideally, 7 out of the 8 
assumptions would be relaxed, making a fit with 15 unknown parameters. This is 
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the full 2-amplitude method. However, if a fit with 15 unknowns is not possible, one 
can reduce the number of unknowns by adding assumptions and incurring a theory 
error. The point is that one has a choice in the type of analysis that is performed. 
In order to maximully reduce the effect of PP, one should include as many free 
parameters in the fit as possible, i.e. make the fewest possible assumptions. 

Now, above we asked the hypothetical question: suppose a future measurement 
finds (f) c s cs = 0.1 ± 0.01 rad (error is statistical only). This disagrees with the SM 
prediction (2/3 s = 0.03636±0. 00170 rad) and therefore suggests NP. Can we conclude 
that NP is indeed present when the theoretical error on (f) c s cs is taken into account? 
We noted that it is estimated theoretically that |A 2/ AAi| = O(10 -3 ). Suppose it is a 
little larger, say 0.005. And suppose that, in the 1-amplitude method, the theoretical 
error is even larger, say ~ 0.04, because many assumptions are made. Thus, taking 
all errors into account, one cannot conclude that NP is present. On the other hand, 
in the 2-amplitude method, the theoretical error is still sufficiently small (~ 0.005) 
that the presence of NP can be confirmed. We therefore see that there are real 
advantages in the search for NP to employing the 2-amplitude method. 

If it is established that the measured value of <fif s differs sufficiently from the SM 
prediction that one can conclude that NP is present, this NP must be in the mixing. 
One can also independently test whether there is NP in the decay. One uses the 
information about the magnitudes and relative strong phases of the Ah and Ah (14 
experimental results + 1 theoretical input). In order to do so, an assumption must 
be made regarding the form of the decay amplitude. For example, suppose that it 
is assumed that NP in the decay is not present, and that the amplitude has the 
form of Eq. ([7]). Now there are 15 unknowns - the magnitudes and relative phases 
of Ai h and A 2t h (the value of the weak phase 7 can be taken from independent 
measurements). The values of these unknowns can be determined from the results 
on Ah and Ah- One can then test the theoretical prediction that |yl 2 /j4i | = O(10~ 3 ). 
If there is a strong discrepancy with this prediction, this might suggest the presence 
of NP in the decay. 

Suppose one assumes that the SM PP term is small, but that there is NP in the 
decay. In this case, the decay amplitude takes the form 

A Jm = Al + e W*pA NP , (23) 

where 0tvp is the (helicity-independent) NP weak phase. But now there are 16 
unknowns - the magnitudes and relative phases of A\ t h and A^p,hi and 4>np- These 
cannot all be determined from the 15 magnitudes and relative strong phases of 
Ah and Ah- However, it has been argued that, to a good approximation, the NP 
strong phases are all negligible [2]. In this case, the A^p^ helicity amplitudes all 
have the same phase (0), and the total number of unknown parameters is 13 - the 
magnitudes of A x h and A NPh (8), the relative strong phases (4), and 4>np- These 
can be determined from the 15 results, so that we can measure the values of the NP 
parameters (given the initial assumption). 
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Above, we have concentrated on B® — > J/ip<fr, but this 2-amplitude method can 
be applied to other B® — > V1V2 decays. For example, consider B° s — > D* + D*~. Its 
amplitude is given by 

A n* a +nr = A w r + {x (s )p , + A ( s)jp , + x u K) + 2 x{ s )pl c w 

= + p' ct - \p' E c w ) + \ { :\PL t - Ip'Fw) 

= Ai + e^A 2 . (24) 

Above, T' and are the colour-allowed tree and colour-suppressed electroweak 
penguin diagrams, respectively. As before, the A 2 term represents the PP, and a 
rough estimate is | ^4 2 / I = O(10 -2 ). If the PP is neglected, the indirect CPA in 
B° s — > D* + D*~ cleanly measures <j)^ cs . However, there is a theoretical error due to 
its neglect. 

This error can be reduced as described above. One does a time-dependent angu- 
lar analysis of the decay, but allows for the possibility of more than one contribution 
to the amplitude, i.e. Ah 7^ Ah- The one difference is that, in B° s — > D* + D*~ , there 
is no 5" helicity. Thus, the index h takes only three values: 0, ||,_L. Consequently, 
there are only six /ifc(t)'s in Eq. (Q. But apart from this, the analysis is unchanged 
- (p c s cs can be measured even in the presence of PP, with a minimal theoretical error. 

Another (slightly different) decay to which this method can be applied is B° s — > 
K*°K*°. This is a pure penguin decay, whose amplitude is given by 

Ar**°=^pi+^K+WPL 

= ^ s) p'ct + ^pu , (25) 

where P' lt = P' u - P[, P' ct = P' c - P[. It is expected that \P' ut \ and \P' ct \ are of 
similar size. Also, in contrast to B° s — > J/ip<ft, these penguin amplitudes are not 
OZI-suppressed. Now, | = 0(A 4 ) and |A^| = 0(A 2 ), where A = 0.23 is the sine 
of the Cabibbo angle. This suggests that the A« ' ] P' uc term can perhaps be neglected. 
However, for consistency, one must neglect allO(X 4 ) terms. One of these is Im(A^), 
so that A^ is real. But since j3 s is related to Im(A| s ' ) )/Re(A^' ) ), it too vanishes in this 
limit. The net effect is that, if one neg lects Ai s) P^ c (the ' 'penguin pollution"), along 
with the other 0(A 4 ) terms, the indirect CPA in B° s — > K*°K*° vanishes because /3 S 
must also be neglected. 

That is, if one wishes to analyze the measurement of B Q S -B Q S mixing in B° s — > 
K*°K*°, the PP term must be kept in the amplitude. Of course, there is now 
a theoretical error in this result, related to the nonzero PP term. In order to 
estimate this error, it has been suggested to measure the corresponding PP term 
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in B° d -»■ K*°K*° and relate it to B° a -»■ ir ^* using flavour SU(3) PSJ. (This 
method also applies to B ds — > K°K°.) As always, there is an uncertainty due to 
the unknown value of SU(3) breaking. 

Once again, one can reduce the problems through a time- dependent angular 
analysis of 5° — > K*°K*°. The method is identical to that in — > J/ip(p, but 
with only three helicities (0, ||,-L), or six /ifc(t)'s. Since both contributions to the 
amplitude [Eq. ( 125]) ] are retained, one can extract 2f3 s with a small theoretical error. 

As is clear from the above, the 2-amplitude method for extracting 2/3 s with a 
minimal theoretical error in the presence of penguin pollution can be applied to all 
B Q S /B° S — > V1V2 decays. (In the case of B° s — > <jxf), in which each <p is detected via 
— > K + K~, it may be necessary to add a fifth helicity. This corresponds to the case 
where both K + K~ pairs are observed in S'-wave.) It can even be used for B Q S /B® 
decays that are governed by a b — > d transition. An example is B® — > J/ipK*°/B® — > 
J/ipK*° [6]. Here, in order that both B° s and B° s can decay to the same final state, 
so as to have an indirect CPA, one must consider the flavour-nonspecific decay of 
the R*°/K*° to K s n . 

To summarize: in the standard model (SM), the weak phase of B° s -B° s mixing, 
2/3 s , is expected to be very small. For this reason, its measurement is of great 
interest, as a disagreement with the SM prediction (2/3 s ~ 0.036 rad) will reveal 
the presence of new physics (NP). The most common decay used for measuring 2f3 s 
is B® — > J/xfxf). Since this decay involves two final-state vector mesons, a time- 
dependent angular analysis must be used to disentangle the CP + and — final 
states. The LHCb Collaboration has determined <j) c s cs (the measured value of 2/3 s ) 
to be small, in agreement with the SM, but still with significant statistical and 
systematic errors. 

In fact, there is also a theoretical error associated with this measurement. B® — > 
J/il><p is dominated by the colour-suppressed tree diagram, but there are other contri- 
butions due to gluonic and electroweak penguin diagrams. As these have a different 
weak phase from the tree, their inclusion in the amplitude leads to a theoretical 
error in the extraction of 2/3 s from the data. They are therefore often referred to as 
"penguin pollution" (PP). In the SM, it is estimated that the ratio of the sizes of 
the two contributions is O(10 -3 ), so that the PP is negligible. Still, there is some 
uncertainty as to its exact size. This is important for the following reason. Suppose 
a future measurement finds <p c s cs to be small, but in disagreement with the SM, even 
including the statistical and systematic errors. Unless the theoretical error is under 
good control, one cannot conclude that NP is present. 

As J/ip and each have spin 1, the B® — > J/i/j4> amplitude is in fact 3 decays: 
the two spins are either parallel (helicity 0) or transverse (two possibilities: helicity 
parallel (||) or perpendicular (_L)). In addition, since (ft — > K + K~, one must also 
allow for the possibility that the observed K + K~ pair has relative angular momen- 
tum / = (S- wave - helicity S). The full B® — > J/ip(f> amplitude therefore involves 8 
amplitudes - Ah and Ah (h = 0, ||, _L, S). When PP in the amplitude is entirely ne- 
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glected (the "1-amplitude method"), one has Ah = Ah. But this really corresponds 
to 8 theoretical assumptions - both the magnitudes and the phases of Ah and Ah are 
taken to be equal. While the theoretical error associated with each assumption may 
be small, the net theoretical error due to the combination of multiple assumptions 
is in fact quite a bit larger. 

In this paper we propose a modification of the above angular analysis in which 
the PP is included (the "2-amplitude method"). The amplitude is written in terms 
of all 8 Ah and A^. We show there are enough time- dependent measurements to 
extract the magnitudes of Ah and Ah and 6 of the 7 relative phases. For the final 
relative phase, theoretical input is required. However, this is not difficult to find. For 
example, we can assume that the Aq-Aq relative phase is zero. With this assumption, 
(pf s can be extracted. The key point here is that the 2-amplitude method only 
involves one assumption, and so the theoretical error is quite a bit smaller than 
that of the 1-amplitude method, which involves 8 assumptions. The 2-amplitude 
method is therefore more promising in the search for NP through the measurement 
of B®-B® mixing. Indeed, if cj) c s cs differs from the SM prediction, this points to NP 
in the mixing, and there is also enough information to test for NP in the decay. 

This same method can be applied to other B° s /B° s — > ViV 2 decays (sometimes 
with slight modification): e.g. B° s D* + D*~, B° s K*°K_*°, B] 00, etc. It 
can even be used for B° s /B° s decays that are governed by a b — > d transitions, e.g. 
B° s J/4>K*°/B° S J/tpK* , in which the K*°/K*° both decay to K s n°. 
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Note added: After this paper was completed, we were made aware of a paper [16] 
that has some overlap with our work. However, most of what is presented in the 
present paper is new. Ref. [IB] concentrates specifically on the case where the tree 
(^4i) and penguin (A 2 ) amplitudes are separated, as in Eq. ([?]). There it is noted 
that, if one takes the mixing phase as known, the angular analysis gives enough 
information to extract 7, given one piece of theoretical input. This is similar to one 
of the points made in this paper: if 7 is assumed to be known, and given one piece 
of theoretical input, the angular analysis can be used to extract not only the B Q S -B® 
mixing phase, but also and A 2y h- However, this is only a minor observation - it 
is discussed in the paragraph preceding Eq. (123]) . The bulk of our analysis concerns 
the extraction of <p c s cs when penguin pollution is present, but not separating the Ai 
and A 2 pieces. That is, we consider the full helicity amplitudes A h and A h [e.g. 
Eq. (1141) ] - the weak phase 7 is irrelevant. Furthermore, we specifically address how 
LHCb's analysis of its data can be modified to incorporate the 2-amplitude method. 
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